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ABSTRACT 



Confexf. The influence of magnetic fields on stellar core collapse and explosion is not well explored. It depends on the possibility to 
amplify the initial pre-coUapse fields. In the absence of rotation this can happen by compression, convection, the standing accretion 
shock instability, and the accumulation and growth of Alfven waves in the accretion flows. 

A\ms. We investigate such amplification mechanisms of the magnetic field during the collapse and post-bounce evolution of the core 
of a non-rotating 15 Mq star with varied initial field strengths, taking into account the microphysical equation of state and neutrino 
physics that play a crucial role in supernova cores. 

Methods. We perform simulations of ideal magnetohydrodynamics with neutrino transport in axisymmetry. The transport of electron 
neutrinos and antineutrinos is treated with a new scheme that solves the energy-dependent set of radiation energy and momentum 
equations in two dimensions by using an analytic closure relation. 

Results. The magnetic field undergoes kinematic amplification by turbulent flows. We also find indications for amplification by 
interacting waves travelling upwards and downwards inside accretion streams. The fields can reach up to equipartition with the 
velocity field. Very high magnetic field strengths require very strong pre-coUapse fields and are able to shape the post-bounce flow, 
leading to a pattern dominated by low-order multipoles. Such models are closest to a successful explosion. 

Conclusions. Magnetic fields can build up to interesting strengths even in non-rotating collapsing stellar cores. Starting with fields in 
the pre-coUapse core as predicted by present stellar evolution models (lO'-lO'" G), typical neutron star fields of lO'^-lO'^ G emerge, 
whereas progenitor fields of some 10"-10'- G lead to fields of magnetar strength of about lO'^'-lO'^ G. Only in the latter case the 
magnetic fields have significant dynamical effects on the flows in the supernova core and may have an influence on the explosion 
mechanism. However, in none of our 2D simulations, we find an explosion until 500 milliseconds post-bounce. 

Key words. Magnetohydrodynamics (MHD) - Supernovae: general - Stars: magnetic fields - Stars: magnetars 
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1. Introduction 

Most scenarios for the explosion mechanism of core-collapse su- 
pernovae (SNe) involve a combination of energy deposition in 
the matter surrounding the nascent proto-neutron star (PNS) and 
multi-dimensional hydrodynamic flows. Examples for means of 
energy transfer to the SN ejecta are the prompt bounce shock, 
neutrinos, magnetic fields and acoustic waves. In the neutrino- 
heating mechanism neutrinos tap the gravitational potential en- 
ergy released during collapse and deposit a part of it behind 
the stalled shock. This process is enhanced and thus supported 
by non-radial fluid flows triggered by hydrodynamic instabili- 
ties like convecti on and the standing accretion sh ock instability 
(S ASI; iBlondin et'al]l2003t rpoghzzdlMll l2002h 

Though it is safe to assume that the exploding star will pos- 
sess a magnetic field of some (unknown) strength and topology, 
its possible influence on the explosion is less clear. The main 
reason for the small number of conclusive investigations into 
this topic is that important effects are expected to occur only 
when the magnetic field is roughly in (energetic) equipartition 
with the gas flow, a condition corresponding to extremely strong 
fields similar to those observed in magnetars. This makes full 
MHD simulations including a treatment of the important neu- 
trino physics in the SN core indispensable. 

Stellar evolution calculations, predicting only weak pre- 
collapse magnetic fields, render the prospects for magneti- 



cally affected explosions very much dependent on the amount 
of field amplification happening during and after collapse. 
Rapid rotation may amplify a weak seed field quite gener- 
ically to dynamically relevant values, e.g., by winding up 
a poloidal field (a process linear in time) or, exponen- 
tially in time , by the magneto-rotational instability (MRI, 
lAkivama et alj 120031) . Magneto-rotational explosions, theoret- 
ically discussed by iMeier et al. (1976), have been studied 



in var i ous approx i mations, e.g., by iBi snovatvi-Kog an et al 



(!l976l);ISvmbalistvl (!l984l);lAkivama et al. (2003); Kotake et aL 
(2004); Obergaulingeretal. (2006); Cerda-Duran et al. (200^ 
Obergaulinge r et al.i (120091) ; recent simulations employi ng de- 
tailed microphysics have been performed by Burro ws et al.l 
(l2007h : IScheidegger et alJ (12008). 

However, according to stellar-evolution models, the major- 
ity of progenitors is expected to rotate slowly. Although many 
stars on the upper main sequence show a surface rotation pe- 
riod close to the critical value for mass shedding, they will most 
likely lose most of their angular momentum during their subse- 
quent ev olution, e.g., by strong stellar winds or magnetic brak- 
ing (Heg er et"al] |2b05; Mevnet et al. 2011). Non-rotating cores, 
albeit lacking the above-mentioned efficient channels for ampli- 
fication, may experience field growth by several effects: 

1. After bounce, convection develops in the PNS and the 
surrounding hot-bubble region. Breaking down into three- 
dimensional turbulence, convection may provide the a effect 
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2. 



responsible for a small-scale dynamo amplifying the field 
on length scales comparable to the scale of turbulent forc- 
ings i.e., the size of convective eddies (Thompson & Duncan 
I1993I) . However, the generation of a large-scale field, e.g., 
on the scale of the PNS, probably requires a non-vanishing 
kinetic helicity of the turbulent flow, which can most nat- 
urally be accounted for by (differential) rotation (see, e.g., 
Brandenbu rg & Subramanian 2005). 
Endeve et al.i ( 1201 Ol) have demon strated that the stan ding- 



Table 1. Variables and symbols used in the MHD equations and 
the neutrino transport formulation presented in this article. 



accretion-shock instability (S ASI iBlondin et al.l 12003 ^. ex- 
cited by a feedback cycle of acoustic and advective perturba- 
tions travelling between the shock wave and the deceleration 
region above the PNS, has the potential to amplify the mag- 
netic field by up to four orders of magnitude. Dynamically 
relevant field strengths (in which case the field reaches at 
least 10% of the equipartition strength) can be reached only 
when the pre-collapse field in the stellar core is sufficiently 
strong. 

3. Non-radial fluid motions triggered by these instabilities can 
excite perturbations of the magnetic field propagating as 
Alfven waves along the field lines. The outward propagation 
of Alfven waves excited close to the PNS has to compete 
with the accretion of gas towards the centre. Assuming that 
the accretion flow decelerates in this region continuously, 
iGuilet et"aD (l2010l) argue that there must be an Alfven point 
at which the Alfven speed equals the accretion velocity, and 
the propagation of the wave (measured in the lab frame) 
comes to a rest. They show that Alfven waves are amplified 
exponentially at such a stagnation point. At the conditions 
of a supernova core, the amplification should be most effi- 
cient for a magnetic field strength of a few 10'^ G and could 
yield final fields of the order of 10^^ G. Dissipation of the 
wave energy can increase the entropy of the gas, thus modi- 
fying the dynamics in the accretion region. For even stronger 
fields, the Alfven point can be close to the shock wave. In 
this case, I Suzuki et al.l (i2008) find that the explosion can 
be driven solely by the energy deposited by the dissipation 
of Alfven waves. This process transmitting energy from the 
(convectively active) PNS to the much less dense surround- 
ing medium bears a strong similarity to the proposed mecha- 
nism for heating the solar corona by Alfven waves emerging 
from the solar surface. 

Previous simulations of magnetohydrodynamic stellar core 
collapse have used a wide variety of methods to treat the 
effects of neutrinos: in the simplest models, t hey were ei- 
ther ignored (e.g., lOber^auhnger et all 120061 : iMikami et alJ 
I2OO8 ) or trea ted by simple parameteri zations or local source 
terms (e.g., ICerda-Duran et all l2008h : more complex ap- 
proa ches include trapping/leakage schemes (e.g., iKotake etaP 
l2004i). mu lti-dimensional, energy-dependent flux-limited diffu- 
sion (|Dessart et al. 20()^ or the isotropic-diffusion source ap- 
proximation which dis criminates between trapped and free- 
streammg components dLiebendorfer et alj 12009,) . In our ap- 
proach to radiation-magnetohydrodynamics (RMHD), we em- 
ploy a new multi-dimensional and energy-dependent scheme 
for the neutrino transport in supernova cores, namely a 
two-moment solver for the neutrino energy (lepton number) 
and mo mentum equations with an analytic closure relation 
(Cernoh orskv & van Weer3ll992l: iPons etalj|2000l: lAudit et alj 
[2002 ). Presently it includes only electron neutrinos and antineu- 
trinos and a relatively simple treatment only of the most relevant 
interaction rates of these neutrinos with neutrons, protons, and 
nuclei. Nevertheless, in the context of the questions focussed on 



Symbol 


Variable 


P 


gas density 




electron fraction 


V 


gas velocity 


b 


magnetic field 


E 


internal energy density 




total energy density 


P 


gas pressure 


P* 


total pressure 




gravitational potential 


rO 



neutrino-matter energy exchange 


rO 
'J v;n 


neutrino-matter lepton-number exchange 


•J V 


IICLILIIIIO liiaLLCi illOlliClilLllii CACilallldC 


Ey 


neutrino energy density 




neutrino number density 


K 


neutrino energy flux 


v 


neutrino number flux 


Pv 


scalar Eddington factor 


p-j 


neutrino energy pressure tensor 




neutrino number pressure tensor 


/ 


0* -moment opacity 




P'-moment opacity 



in this paper, the scheme provides a reasonably good representa- 
tion of the neutrino effects that play a role during core collapse, 
bounce, shock propagation and in the accretion layer behind the 
stalled supernova shock. Moreover, it is a computationally effi- 
cient treatment of the neutrino transport, and two-dimensional 
simulations up to several hundred milliseconds after bounce are 
easily possible. 

It is the goal of this work to investigate the relevance of field 
amplification mechanisms like those in points 1.-3. concerning 
their importance for the evolution and potential revival of the 
stalled shock by core-collapse simulations including a reason- 
ably good treatment of the relevant microphysics. In particular, 
we will study the field growth connected to convective and SASI 
activity in the post-shock layer and the role of energy transport 
and dissipation by Alfven waves. To this end, we perform MHD 
simulations of collapse and post-bounce evolution of the core 
of a star of 15 solar masses (Wooslev et al. 2002). Varying the 
strength of the initial field, assumed to be purely poloidal, we 
can determine different sites and mechanisms of field amplifica- 
tion and identify the back-reaction of the field onto the flow. 

This article is organised as follows: Sect.|2] describes our 
physical model and the numerical methods; Sect. [3] introduces 
the initial conditions; Sect.|4] presents the results of our simu- 
lations; Sect.|5] gives a summary of the study and draws some 
conclusions; appendix|A]gives a brief overview of our treatment 
of the neutrino transport. 



2. Physical model and numerical methods 

2.1. Magnetohydrodynamics 

We assume that the evolution of the gas and the magnetic field 
is described by the equations of Newtonian ideal magnetohydro- 
dynamics (MHD), 



d,p + V ■ (pv) = 0, 
5,(pye) + V ■ (pFeV) = Sl„ 
d,{pv') + Vj [P^6'j + pv' v> - b'bJ) = pV'<D + S\„ 



(1) 

(2) 

(3) 
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d,e^ + V ■ ((e* + n)v - (V ■ b)b) = pv ■ V(D + 5[ 
- V X (v X ft) = 0, 



(4) 
(5) 



describing the conservation of mass, electron fraction, gas mo- 
mentum, total energy of the matter, and magnetic flux, respec- 
tively. In addition to these evolutionary equations, the magnetic 
field has to fulfil the divergence constraint. 



Table 2. Average magnetic field strengths and normalised mag- 
netic energies for our computed models with magnetic field (ex- 
cluding model slS-BQ). See text for details. 



V ■ ft = 0. 



(6) 



The symbols used in this system have standard meaning; they are 
listed in Tab.[T] The total energy density and pressure are defined 
as the sum of the contributions of the gas and the magnetic field, 
e* = e -H 1 /2 pv^ -H 1 /2 ft^ and = P + 1 /2 ft^ respectively. 

The MHD part of th e system is closed by the equation of 
state of IShen etal ] (fT998h . relating gas density, internal-energy 
density, and electron fraction to temperature, T, gas pressure P, 
entropy, S , the composition of four species (neutrons, protons, a 
particles, and a representative heavy nucleus), and the chemical 
potentials of neutrons, protons, and electrons. 

To include some key features of general relativity into our 
Newtonian model, we use the modified gravitational potential of 
iMarek et al.l (l2006l) . Following the lines proposed by these au- 
thors, we include the gas and neutrino energy density and pres- 
sure as well as an approximate Lorentz factor of the flow in the 
source term of the Poisson equation for <t). Overall, this leads to 
a somewhat deeper potential well when the core approaches nu- 
clear matter density and serves as a good approximation to the 
stronger relativistic gravity of a PNS. 

Our code is based on an implementation of the MHD 
equations on a fixed Eulerian grid in flie constrained- 
transport discretisation (lEvans & Hawlevlll988l) . It uses high- 
resolution shock-capturing (HRSC) methods with approx- 
imate Riemann so l vers (Lax-Friedrichs, HLL, or HLLD 
(iMivoshi & Kusand l2005h ) n the multi-stage framework 
dTitarev & Torol 120051) . We reconstruct the conserved vari- 
ables using high-order monotonicity-preserving (MP) methods 
dSure sh & Huynh 1997). Time stepping is done by explicit 
Runge-Kutta schemes of 2"'*, 3"*, or 4* order. 



2.2. Neutrino transport and neutrino-matter interactions 

We treat the transport of neutrinos by solving the energy- 
dependent, multi-dimensional moment equations for neutrino 
energy, Ey, and neutrino momentum, f J,, 



d,Ey + V ■ (EyV + Fy) - OjV jV kO ,,P'y ^ S^, 

drF'y + Vjify +p'j) + FlVjVj = S'y'\ 



(7) 
(8) 



Presently only electron neutrinos, Ve, and antineutrinos, Ve are 
considered. The source terms, 5^', include the interaction of 
neutrinos with the gas by emission, absorption, and scattering 
reactions; we use a reduced set of processes, viz. 

1 . emission and absorption of electron neutrinos by neutrons, 

2. emission and absorption of electron anti-neutrinos by pro- 
tons, 

3. elastic scattering of all neutrinos off nucleons, 

4. emission and absorption of electron neutrinos by heavy nu- 
clei, 

5. coherent elastic scattering of all neutrinos off heavy nuclei. 

Furthermore, fluid-velocity dependent effects such as the P dV 
work associated with diverging flows appear in the equations 
as well. To close the system of moment equations, we have to 



model 


slS-BlQ 


slS-Bll 


slS-Bll . 5 


S15-B12 


bo 
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1.4 X 10'" 
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7.8 X 10"' 
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700 
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2.1 X 10"* 


4.7 X 10'" 




1.1 X 10-' 


2.9 X 10-5 


9.8 X 10-5 


4.8 X 10-" 


0" 

"cnv 


0.00041 


0.055 


0.16 


1.6 


Z^stb [G] 


3.8 X 10'- 


3.0 X 10'' 


6.3 X 10'' 


1.4 X 10'" 




1.1 X 10"* 


7.8 X 10-5 


3.3 X 10-" 


1.4 X 10-' 




0.0020 


0.16 


0.80 


3.7 


Z^gain [G] 


3.6 X 10" 


6.0 X 10'" 


1.0 X 10' ' 


1.2 X 10' ' 




2.2 X 10"* 


5.7 X 10-» 


2.7 X lO-'' 


1.0 X 10-2 


/ 


0.000089 


0.016 


0.071 


0.26 



specify the tensor of the second mom ent, Pj/, for whic h we use 
a simple analytic closure relation from iMinerbol d 1 97 8l) . Using a 
tensorial generalisation of a one-dimensional Eddington factor, 
our method is generically multi-dimensional. A more detailed 
discussion can be found in Appendix lAl 

We note that our set of equations is a generalisation of the 
usual diffusion ansatz, in which the system of moment equations 
is truncated at the level of the energy equation and closed by 
expressing the flux in terms of the gradient of the energy density, 
mostly connecting diffusion and free-streaming limits by a flux 
limiter. Retaining the first two moments, it leads to a hyperbolic 
system, which can be solved by common methods such as high- 
resolution shock-capturing methods. More details on this can be 
found in Appendix lAl 



3. Models and initial conditions 

Interested in the contribution of magnetic fields to the dynamics 
of supernovae, we simulate the evolution of a star o f 15.0 so- 
lar masses with solar metalicity dWooslev et al.ll2002l) . We map 
the pre-collapse model to a grid of - 320 radial zones dis- 
tributed logarithmically over a radial domain of [r^in = 0; rmax ~ 
2800 km] and ng = 64 lateral zones distributed uniformly be- 
tween the north and south pole. The resolution at the grid centre 
is (5r)cti- - 400 m. Ten energy bins between emin =0.1 MeV and 
Emax - 280 MeV are used for the neutrinos. 

The topology of the magnetic field at the onset of collapse 
is highly uncertain. On the main sequence, field amplification 
by, e.g., gradual contraction of the star or convection competes 
with the loss of magnetic energy in stellar winds and in work 
the magnetic field does by exerting torques on the stellar mat- 
ter. In the absence of rotation, stars lack an important ingredient 
of most large-scale dynamos. In such a case, they may be dom- 
inated by a small-scale turbulent field rather than a large-scale 
field, e.g., a lower-order multipole. Nevertheless, we assume a 
simple initial field, viz. a modified dipole (Fig.lTJ. While this 
is probably not a very good approximation to real stellar cores, 
it represents an optimal configuration for the Alfven-wave am- 
plification mechanism we are interested in because of the large 
coherent radial component favouring the radial propagation of 
Alfven waves. Furthermore, we observe that the field is replaced 
by a more complex small-scale field in the regions of the core af- 
fected by hydrodynamic instabilities such as convection and the 
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log Ibl [G] 

7.92 8.63 9.33 10.04 




X [1000 km] 

Fig. 1. Initial (pre-collapse) field configuration of model 
slS-BlO: field strength (colour scale) and field lines. For sim- 
plicity, only one quadrant is shown. 

SASI. Hence, we deem the influence of such an artificial choice 

for the initial field on our results only a min or one. 

The field is the same as the one used bv lSuwa et al ] (l2007h . 
defined by a vector potential of the form 

A^ = ^^^r^in6. (9) 

We set the normalisation radius ro denoting the location of the 
dipole to ro - 1000 km and vary the parameter /jq setting the 
field at the centre of the core between 10 '° and 10'^ G, i.e ., in 
a range somewhat above that found by (iHeger et al.ll2005l) for 
rotating progenitors. See Fig. [T] for a visualisation of the initial 
field of model sl5-BlQ. 

Tab.|2]list important parameters of our models: apart from 
the initial field strength, b^, we present time averages (over t e 
[280 ms, 300ms]) of the mean field strength in the PNS at den- 
sities above lO'"* gcm"^, in the convection zone inside the PNS 
(^cnv), in the stable layer surrounding the PNS (/^stb), and in the 
gain region where neutrinos deposit energy behind the SN shock 
(^gain)- Moreover, time-averaged (over t e [280 ms, 300 ms]) ra- 
tios of magnetic energy to internal energy and to kinetic energy 
in the same regions, f? and fi^, respectively, are listed. We will 
define these regions more quantitatively in Sect. l4.1l 

4. Results 

4.1. A non-magnetised reference model 

The model with bo - 0, model sl5-B8, serves as a reference 
case to which we can compare the magnetised models; further- 
more, it allows us to compare our results with previous results 
for similar cores. The evolution of the mean entropy profile and 
of the root-mean-squared lateral velocity of the model as func- 
tion of time is shown in Fig.|2] 



Core bounce occurs after 168 milliseconds of collapse; the 
shock forms at a mass coordinate of M^h ~ 0.45 Mq and starts 
propagating outwards. A strong negative gradient of the elec- 
tron fraction triggers prompt PNS convection in a layer between 
an inner radius of about 15 and an outer radius that is initially 
at about 50 km and decreases to about 40 km at the end of the 
simulation. In Fig.|2] we show the upper boundary of the PNS 
convection zone by a solid black line. Convective velocities in- 
side the PNS reach maximum values of » 5 x 10** cms"' and 
angle-averaged rm.s. values of x 2.5 x 10^ cms"'. 

After roughly 50 ms post-bounce, the shock st arts to oscillate 
as the standing accretion shock instability (SASI Blondi n et al.l 
I2OO3.) develops in the post-shock region. Perturbations created 
at the shock wave are advected with the flow towards the PNS, 
leading to the parallel alignment of the patterns of lateral ve- 
locity and of the (grey) mass-shell lines lines in Fig.|2] From 
these perturbations, the SASI develops in the post-shock region. 
Typical transverse velocities associated with the SASI reach val- 
ues about 3 times as large as the convective velocities deeper 
inside the PNS. With increasing time, the SASI begins, starting 
in the immediate post-shock region, to affect a larger fraction of 
the post-bounce volume, until at f > 100 ms the PNS convection 
zone is separated from the SASI region by a thin stable layer 
located roughly between densities of lO" and lO'^ gcm"^. The 
border between the stable layer and the SASI region is shown in 
Fig.|2]by a solid black line. 

For later reference, we list the definitions we use for the four 
distinct regions of the core and their respective boundaries: 

1. The inner hydrodynamically stable PNS extends from the 
origin to roughly the radius where the density of the gas 
drops below lO''* gcm"^. 

2. This inner core is surrounded by the PNS convection zone, 
characterised by a negative gradient of the electron fraction. 
The outer boundary of this layer is associated with the loca- 
tion of the minimum of the profile. 

3. A stable layer surrounds the PNS convection zone and sep- 
arates the latter from the 

4. gain region where hot-bubble convection and the SASI op- 
erate. The boundary between the stable and the gain layer 
is approximately given by the gain radius that defines the 
transition from neutrino cooling below to neutrino heating 
above. 

Despite growing SASI activity and increasing entropy in 
the post-shock region and despite a slowly growing radius of 
the stalled shock, no explosion takes place until 550 ms after 
bounce. This is in agreement with simulations for 15Mfr, stars 
(b ut different p rogenitor models) by .Marek & Janka ( 2009.) and 
iBurrows et al.l(l2006l) . 

4.2. Models with weak initial fields 

Models S15-B1Q and sl5-Bll (^o = 10'°'" G, respectively) 
exhibit very similar dynamics compared both to each other and 
to the non-magnetised reference model. The structure of the core 
is during several hundred milliseconds of post-bounce evolution 
the same as that of model sl5-B(S: a convective layer inside the 
proto-neutron star at p > 10'^ gcm"^ and a SASI region, sep- 
arated by a thin stable deceleration layer; hence, the profiles of 
entropy and vq are in qualitative agreement with the ones shown 
in Fig.m 

We discuss the evolution of the magnetic field of model 
S15-B10 in the following. A space-time diagram of the model 
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t[nii] 

□ im ym -ma $ud 




a im '^3lnf ^an son 

trmii 

Fig. 2. Summary of the evolution of model sl5-B0. The upper and lower panels show angular averages of the entropy per baryon 
and the of root-mean-square of the lateral velocity, respectively (colour coding) as a function of radius and post-bounce time. The 
grey lines indicate the positions of different mass shells; the intervals between solid and dashed lines are 0.1 and 0.025 solar masses, 
respectively. The thick black lines show (from small to large radii) the outer boundaries of the PNS convection zone and of the 
stable layer beneath the SASI region and the outermost position of the shock wave. The thin dashed black lines depict iso-density 
surfaces of 10" gcnr^,n - 14, 13, .... 



(Fig.O shows angular averages of the current densitjQ j - 
(V X ft)* and the field strength. We use the current density to 
detect spatial variations of the magnetic field. While not all pat- 
terns of j are necessarily caused by Alfven waves, Alfven waves 
are associated with a non-zero current density. Hence, we can try 
to identify these waves in the current distribution. 

Without changing the field geometry, compression amplifies 
the magnetic field during collapse along with the maximum den- 



' For simplicity, we omit the factor ^ in the definition of the current 
density and consider only the 0-component, neglecting the vanishing 
components /' and /. 



sity as Pmax' leading to a gain of a factor of ^ 500 and ^ 25 in the 
maximum field strength and total magnetic energy, respectively. 

Further amplification occurs in the PNS convection zone (see 
Fig.lUi. Within roughly 10 ms after the onset of PNS convection, 
the mean field is amplified by a factor of ~ 3 between radii of 
~ 20 and ~ 60 km. Field growth is driven by the convective over- 
turn motions of the gas (indicated in the deviations of from its 
angular mean value, see right half of Fig.|4]i and leads to a fil- 
amentary field topology (same figure, left half). This strongly 
non-uniform structure is reflected in the large fluctuations of the 
current density visible in the top panel of Fig.[3](first ~ 100 ms 
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IDD ma ^^sfftf; 4m 5m 

trmil 

Fig.3. Evolution of model sl5-B10. The panels show in colour coding the logarithm of the lateral average of the current density 
(top panel; the intensity of the colour encodes the absolute value and the blue and red colours distinguish between positive and 
negative signs) and the logarithm of the mean magnetic field strength (bottom panel). The lines have the same meaning as in Fig.|2] 



at radii between ~ 15 and ~ 60 km). After about 100 ms, the 
amplitude of variations of the current density in the convective 
layer decreases. While j still varies on short temporal and spa- 
tial scales, larger coherent fluctuations are mostly absent during 
this phase: compare the larger red (j'^ < 0) patches in the PNS 
convection zone at f < 100 ms to the convection zone at later 
times standing out as a broad blue (j'^ > 0) band with little con- 
tribution of negative currents. In this phase, the angular average 
of the current density in the PNS convection zone is very regular 
and shows little temporal variation. It is dominated by one large- 
scale pattern, viz. columns of strong magnetic field generated by 
flows converging at the axis due to the constraint of axisymme- 
try (Fig.m. These very persistent features are associated with a 
strong positive current density whose value is dominating what- 



ever fluctuations there may be due to convective cells at the same 
radii. 

As soon as the SASI develops, we find appreciable lateral 
velocities also in the region above the PNS convection zone. An 
early stage of the development of the SASI can be seen akeady in 
the Ye variation close to the shock in Fig.|4| small deformations 
of the shock excite perturbations advected inwards. These per- 
turbations are associated with kinks in the magnetic field lines. 
Geometrically, they take the form of sheets oriented parallel to 
the shock, i.e., roughly at constant radius, containing a magnetic 
field with a strong fl-component. This translates into a non-zero 
derivative drb'', i.e., a non-zero current density. Consequently, 
the perturbations can be identified easily in the top panel of Fig. [3] 
as the alternating blue and red bands starting at the shock and 
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Fig. 4. Snapshot of model sl5-B18 shortly after the onset of 
PNS convection at r < 50 km. The left and right panels show 
the logarithm of the magnetic field strength and the deviation of 
the electron fraction from its angular average. Vectors indicate 
the velocity field. 



>, 




100 



Fig.5. Snapshot of model sl5-B10 ait k 120 ms post-bounce. 
The panels show, apart from the velocity field, the logarithm of 
the magnetic field strength (left) and the current density (right). 



falling towards the deceleration layer Propagating at the veloc- 
ity of fluid accretion, these bands are parallel to the mass shells 
(grey lines). 

When the accretion flow decelerates above the PNS convec- 
tion zone, these sheets pile up at the deceleration layer (cf. Fig.|5] 
at r > 50 km), leading to a local increase of the magnetic energy. 
The bottom panel of Fig. [3] indicates that this process happens 
throughout the entire evolution and magnetic energy continues 
to accumulate. 



10 



10 



r Pan] 



100 



Fig. 6. Average radial profiles of the field strength of model 
sl5-Bl(5 in the post-bounce phase (time averages over t € 
[100; 200] ms (black), t e [200; 300] ms (green), t € 
[300; 400] ms (red) post-bounce). The thin straight lines show 
profiles oc (r < 50km) and oc r^^^^ ( r > 80 km). 



4.2.1. Sub- and superalfvenic regions 

In the following, we discuss the propagation of MHD waves in 
the core. We will not touch upon fast modes and focus only on 
slow and Alfven waves. Whereas Alfven and slow modes prop- 
agating across the field lines are clearly distinct solutions (they 
have, e.g., different polarisation states and different propagation 
speeds), they degenerate when propagating along the field lines. 
In this case, they have identical propagation speeds, viz. the 
Alfven speed. Because we expect the dynamics to be domi- 
nated by the modes propagating fastest, i.e., the ones travelling 
(anti-)parallel to the field lines, we will consider these modes 
only. Because of the degeneracy of slow and Alfven modes in 
such a case, we will refer to these modes as Alfven waves for 
short, postponing an investigation of the possible differences be- 
tween the two classes of modes. It is expected (J. Guilet, pri- 
vate communication) that slow modes behave very similarly to 
Alfven waves. 

An essential requisite for this evolution is that the flow is 
superalfvenic, i.e., the flow velocity is larger than the Alfven ve- 
locity of the magnetic field 



= 10" cms" 



\b\ 



10'2G 



lO'^gcm- ^^'^' 
P 



(10) 



A magnetic field b allows for the propagation of Alfven waves 
in and against the direction of the field at the same propagation 
speed, ca relative to the flow of the gas carrying the magnetic 
field, i.e., with a velocity +Ca = ±bl\b\ ca- The total velocity of 
such a pair of waves in the lab frame, where the gas moves with 
a velocity v, is v ± ca. Hence, an Alfven point develops where the 
component of the Alfven speed directed along the fluid velocity 
has the same absolute value as v. At that point, one of the pair of 
Alfven waves, viz. the one with an Alfven velocity antiparallel 
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Fig. 7. Space-time diagram of the magnetic field strength (top panel) and the logarithm of the Alfven number (bottom panel) of 
model S15-B11.5 along the north pole. 



to V, is trapped , i.e., has vanishing total velocity, while the other 
one passes the Alfven point at a total speed 2|v|. 

To quantify the importance of the magnetic field, we define 
the total, parallel, and perpendicular Alfven numbers. 



Atot — 
All = 

A^ = 



|caI 
V ■ Ca 

•-A 

IvxcaI 



(11) 

(12) 
(13) 



- A|| regulates the total (lab-frame) propagation speed of the 
wave, composed of the (comoving) Alfven velocity and the 
flow speed, +Ca + v; |A||| = 1 indicates an Alfven point where 
the waves stall. 

- Ax is a parameter relevant to the excitation of perturbations 
orthogonal to the magnetic field, comparing the velocity dis- 
placing magnetic field lines with the tension of the field lines. 
Alfven waves rely on magnetic tension as restoring force; 
hence, their excitation occurs for < 1; otherwise, the 
field can be distorted (e.g., amplified in a stretch-twist-fold 
dynamo), but the distortions are beyond the linear regime of 
Alfven waves. 



While Atot compares kinetic and magnetic energies directly. Ay 
and Aj^ measure two effects important to Alfven waves propa- 
gating along the field: 



The gas flows triggered by PNS convection as well as by 
the SASI are highly superalfvenic. Though there is, in princi- 
ple, a small stable region between the convection and the SASI 



M. Obergaulinger and H.-Th. Janka: Magnetic field amplification in collapsing, non-rotating stellar cores 



9 



regions, where the accretion is decelerated and only small lat- 
eral velocities are present, the flow in this region remains super- 
alfvenic because of overshooting from the PNS convection and 
the SASI regions. Travelling along a field line connecting the 
convection zone with the SASI region, one observes a drop of 
all three Alfven numbers in the deceleration region, but they still 
remain considerably above unity. Hence, there is neither any sig- 
nificant excitation of Alfven waves nor is t here an Alfven point 
where amplification as proposed by Gui let et alj ((2010) could 
happen. Magnetic perturbations are unable to travel upwards 
against the accretion flow; furthermore, the relative weakness of 
the field prohibits dynamic back-reaction. Thus, this is a purely 
kinematic effect. 

The combination of these effects leads to the time evolu- 
tion of the average field profile shown in Fig.|6l Across the PNS 
convection zone, the average field strength drops by an order of 
magnitude. At its upper boundary, the pile-up of field accreted 
from the outer layers leads to a strong local increase of the field 
strength, while in the SASI region, the radial decrease of the field 
continues following a slightly flatter law (b cc r^^/^) than in the 
PNS convection zone (b cc r"^). The steepening of the slope at 
r ~ 200 km marks the minimum radius of the shock wave. 

For model sl5-B10, the flow is superalfvenic almost 
everywhere. Sizeable subalfvenic regions develop only for 
significantly stronger fields. The convective flows of model 
S15-B11.5, with an initial field strength given by bo - 
lO"^ G ~ 3.16 X 10" G, are only moderately superalfvenic. In 
particular, along the pole, where the magnetic field is enhanced 
artificially for symmetry reasons, the Alfven numbers are below 
unity (Fig.|7j- The thin layer separating PNS convection and hot 
bubble consists of both sub- and superalfvenic regions, while the 
accretion flow is mostly superalfvenic. 

In analogy to the model discussed above, perturbations of the 
magnetic field created at the shock wave are advected towards 
the PNS. The field piling up above the PNS, dominated by a 
strong lateral component generated by the complex overturns of 
the gas flow, forms shells around the PNS. Since the radial com- 
ponent is weak, Alfven waves are bound to travel at roughly con- 
stant radius. Consequently, the location of an Alfven point along 
a field line is determined not by the radial flow alone but mostly 
by the lateral velocity. Whereas the former tends to go through 
zero in the deceleration region, the latter can, in principle, have 
a stochastic and chaotic behaviour and dynamical evolution. 

Thus, while there may be Alfven points in this region, their 
location is highly variable; they may even, depending on the 
flow, cease to exist for a certain time when the flow becomes 
superalfvenic or subalfvenic in the entire region. This may have 
adverse eff'ects on the efficiency of the amplification process 
proposed by Guilet et al. (2010). Despite the stochastic nature 
of these variations, we find a systematic trend: the stronger the 
magnetic field is, the farther out stretches the layer of accumu- 
lated magnetic field on average. Since this layer corresponds 
to the hydrodynamically stable region between PNS convection 
and the SASI region, the trend of a growing radius of the region 
where the field accumulates can be seen most easily in the aver- 
age profiles of the lateral velocity shown in Fig.|8] This effect is 
an, albeit small, indication of dynamic back-reaction of the field 
onto the flow. 

4.2.2. Accretion columns 

The accretion columns display much stronger radial components 
of the field and the velocity, while non-radial fields of a consider- 
able strength may be present in the flux sheets surrounding them. 




r [ km] 

Fig. 8. Root-mean-square averages of the fl-component of the 
velocity as a function of radius for models sl5-B10 (black), 
s 1 5 -B 1 1 (red), and sl5-B11.5 (blue). Time averages over sev- 
eral tens of milliseconds in the late post-bounce phase (at t a; 
300 ms) are shown. 



Accretion proceeds mostly through a narrow stream at the poles, 
and mass conservation leads to an increase of the flow speed 
as the gas is accreted through a small solid angle rather than 
uniformly over all angles. Eventually, the flow reaches sound 
speed. Analogously, flux conservation increases the magnetic 
field strength in the accretion funnels. Depending on the ini- 
tial field strength, the Alfven number of the accretion flow and 
/3 = P/Pmug may be of order unity or much greater 

Due to the high Mach number, pressure perturbations of the 
accretion flow may steepen, in some cases yielding discontinu- 
ities. In particular, weak shock waves may mark the lower end 
of the column, providing a very abrupt deceleration of the gas. 
Due to the dissipation of kinetic into internal energy, the entropy 
of the matter above the shock is higher than that of the gas in 
the deceleration and PNS convection layers where S is nearly 
uniform. 

In the presence of a magnetic field, these features are not 
purely hydrodynamic but rather MHD structures associated with 
a bending of the field lines, i.e. a current density: the field, pre- 
dominantly radial in the accretion column, is folded by the rel- 
atively high lateral velocity in the deceleration region and is 
mostly mostly lateral in this layer. This effect is weakest in the 
polar accretion funnels because the axial symmetry restricts the 
development of a non-radial field considerably. In all accretion 
columns, both velocity and magnetic field are dominated by their 
respective radial components; for the polar accretion flow, this 
dominance of the radial components holds even at radii below 
the lower end of the accretion stream, i.e., for r < 70 - 100 km. 
Accretion streams at intermediate latitudes tend to punch into 
gas threaded by strong lateral field lines and moving at consid- 
erable speed along constant radius. 

For a sufficiently strong magnetic field, an Alfven point may 
be located in the accretion stream; for a steep gradient of or 
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even a discontinuity at the lower end of the stream, this is a most 
natural place for the Alfven point. 

Conditions are most favourable for the formation of a stable 
Alfven point in radial polar accretion streams. We show the time 
evolution of the magnetic field and the Alfven number along the 
north pole of model sl5-Bll . 5 in Fig.|7] Initially, accretion is 
roughly spherical, but after f « 80 ms, polar accretion streams 
begin to form; as a consequence, the radial velocity behind the 
shock increases abruptly. This stream remains active for the rest 
of the simulation. It is, however, subject to strong fluctuations 
on time scales of about 10 milliseconds. These oscillations af- 
fect the outer and inner termination radius of the accretion flow, 
the velocity, and the magnetic field. They are caused by lateral 
motions due to convection and the SASl in the surrounding gas 
that lead to an episodically varying cross section of the fun- 
nel. Following the variations of the accretion velocity, the Mach 
number of the flow oscillates around unity. For this model, the 
same holds for the Alfven number 

The accretion flow terminates in a steep gradient hosting an 
Alfven point. This point is first, at f s; 100 ms post-bounce, 
situated at a radius of 80 kilometres, i.e., about 30 kilome- 
tres above the PNS convection layer; later it retreats to lower 
radii as the core contracts. The intermediate layer between the 
termination point of the accretion column and the PNS con- 
vection layer is, though threaded by a weaker magnetic field, 
filled by gas at subalfvenic velocities. Hence, at these radii, 
Alfven waves could propagate from below towards the accre- 
tion flow. Given a typical Alfven velocity in this layer of ca ~ 
2 X 10^ cms~'(feo/10" G), they should traverse the intermedi- 
ate layer within ta ~ 10 ms/(Z7o/10" G), i.e., on the time scale 
of the oscillations of the accretion column. This has two conse- 
quences: first, it complicates the identification of Alfven waves, 
and, second, it may limit the applicability of the scenario of am- 
plification of Alfven waves in this context. 

The fluctuations of the accretion velocity and the position 
(radius, possibly also latitude) of the termination point of the ac- 
cretion column induce perturbations propagating in the layer be- 
tween the accretion column and the PNS convection region, i.e., 
outside a radius of about 50 kilometres and inside a radius that 
recedes gradually from about 80 - 100 kilometres attx 100 ms 
to roughly 60 km at f ~ 300 ms. Partially, these perturbations can 
be entrained in matter flowing upwards and thus sustaining the 
magnetic field in the hot-bubble region. Additionally, overshoot 
from the convective layer below contributes to a certain degree 
to maintaining the turbulence in the deceleration region. 

4.2.3. Field amplification 

Though, in principle, the dynamics in this layer can be expressed 
in terms of basic MHD waves, it is not trivial to identify these 
simple solutions in the highly variable fields, preventing deeper 
insights into the field amplification mechanisms. Therefore, we 
must try to find indirect evidence for the amplification processes 
at work. To this end, we distinguish between kinematic and dy- 
namic amplification, i.e., processes that do not depend on the 
seed magnetic field (at least during an initial phase when it is 
still sufficiently weak) and processes for which the amplification 
is some function of the seed field, respectively. 

Convection and the SASI amplify the field kinematically. 
Hence, they should amplify the field by a given factor indepen- 
dent of the seed field. This is demonstrated by the comparison 
of space-time averages of radial profiles of the magnetic field 
strength of models with different initial fields shown in Fig.|9] 
Corrected for the scaling with the initial field strength, the pro- 
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Fig. 9. Space-time averages of the magnetic field strength of 
models sl5-B10 (black), sl5-Bll (green), and sl5-Bll . 5 (or- 
ange line). We average over times between t = 200 ms and 
t - 300 ms post-bounce. The fields of models sl5-Bll and 
slS-Bll . 5 are scaled according to their initial strength. 
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Fig. 10. Left panel: Temporally averaged (over post-bounce 
times of 120ms < t < 200 ms) radial profiles of the magnetic 
field strength along the north pole for models sl5-B18 (black), 
sIS-Bll (green), sIS-Bll . 5 (orange). The field strength is rep- 
resented by solid lines; the dashed lines show the field strength 
scaled by the ratio of initial field strength. 

files of all three models are remarkably similar. The main differ- 
ences appear in the deceleration region, where neither the radial 
collapse nor convection or the SASl are the main amplification 

mechanisms. 

The Alfven mechanism due to iGuilet et"al] (1201 Oh . on the 
other hand, depends on the presence and the location of 
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Fig. 11. The top panels show profiles of the 
radial velocity (solid lines) along the north 
pole of models sl5-B10 {left) and sl5-Bll. 5 
{right); the dashed lines show the Alfven ve- 
locity, c'A, and its negative, -Ca- The bottom 
panels depict profiles of the magnetic field 
strength of the same models. In all panels, 
black and green lines denote the profiles at 
t X 100 ms and t ~ 200 ms post bounce, re- 
spectively. 
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Alfven points, hence on the initial field strength. Furthermore, 
the Alfven-wave mechanism should be most efficient along the 
polar axis (as discussed above and shown in Fig.|7]i. A compar- 
ison of the polar regions of models with different initial field 
strengths should therefore allow us to determine the importance 
of the Alfven-wave mechanism. 

We compare time averages of the field strength along the 
north pole for different initial fields in Fig.[TO] and profiles for 
two times in Fig.[TT] In the PNS convection zone, the field is 
amplified roughly by the same factor irrespective of b^, but only 
until it reaches kinetic equipartition. This limit is reached by 
model S15-B11.5. The field decreases with radius in the con- 
vection zone and it exhibits a further drop by about an order of 
magnitude at the upper boundary of the convective layer inside 
the PNS. 

After a few tens of milliseconds, the evolution of models 
s 1 5 -B 10 and sl5-B11.5 diverges. The weak-field model shows 
a coherent stable downflow inside the convection region (at radii 
between 20 and 50 kilometres) and a rather smooth velocity pro- 
file in the surrounding layers (at radii above 70 kilometres; see 
the black line in the upper left panel of Fig.fTTTi. Alfven points 
may be located at the border between the PNS convection region 
and the accretion flow. The accretion flow transports the mag- 
netic field towards the convection zone and establishes a tem- 
porally averaged profile varying around |v| oc between the 
convection zone and the shock wave. 

In model sl5-Bll . 5 (for profiles at different times, see the 
right panels of Fig.fTTI). strong perturbations originating in the 
convection zone propagate upwards against the accretion flow 
(note the large positive velocity below r » 50 km). Under the 
influence of perturbations converging from below and above. 



the lower end of the accretion column develops an increasingly 
steep profile until a discontinuity is formed a few tens of kilo- 
metres above the PNS convection zone (visible in the green line 
at r ^ 100 km; the discontinuity located at r » 200 km is the 
stalled supernova shock wave). Magnetic energy is transported 
towards this point from both sides (some of the accreted field can 
be identified in the black curve in the lower right panel above 
r ^ 100 km) and a region of enhanced field correlated with the 
termination of the accretion flow is created, visible right at the 
discontinuity (at r > 90 km in the green line). 

Thus, upflows travelling from the PNS convection zone to- 
wards the polar accretion columns play a role in stopping the 
accretion flow and shaping the radial profile at its lower end. 
We find that the lower end of the equatorial accretion column 
develops steeper profiles, too, but this is caused by a different 
effect: the accretion flow continues to smaller radii until it is ter- 
minated at a radius of r s; 60 km, i.e., it terminates when it 
hits the high-density layer near the outer boundary of the PNS 
convection zone. These discontinuities do typically not coincide 
with an Alfven point and are not associated with an enhancement 
of the magnetic energy. 

Hence, we find that field amplification, while dominated by 
kinematic effects such as convection and the SASI, has also a 
certain contribution from the growth of the field by the super- 
position of perturbations propagating in opposite directions be- 
tween the regions of PNS convection and the SASI. 

4.3. Stronger magnetic fields 

For an initial field of bo = 10'^ G (model sl5-B12), we find, 
much like in the models with weaker fields, a convection layer 
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Fig. 12. Snapshot of model sl5-B12 at f ~ 130 ms post-bounce, 
described roughly, and the top and bottom right halves show the 
logarithm of the parallel and perpendicular Alfven numbers, re- 
spectively. Sub- and superalfvenic regions appear red and blue, 
respectively. 
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Fig. 13. Average (over angle and time, 160 ms < t < 220 ms 
post-bounce) profiles of the magnetic field strength of models 
S15-B12 (black), sl5-Bll (green), and sl5-B10 (orange), cor- 
rected for the scaling of the initial field strength in the case of 
the latter two models. 




inside the PNS and, outside of the PNS, the hot-bubble region 
dominated by the SASl. Due to the much stronger initial field, 
however, the accretion flow onto the PNS is subalfvenic almost 
out to the shock wave. 

The spatial structure of the model at f ^ 130 ms is dis- 
played in Fig. [12] The hot-bubble region is dominated by a super- 
equipartition (w.r.t. the velocity) magnetic field; this is particu- 
larly evident in the large-scale arcs of the magnetic field right be- 
hind the shock wave. In these arcs, which separate gas of a high 
entropy on the inside from low-entropy gas on the outside, the 
magnetic field reaches equipartition with the internal energy of 
the gas. Gas falling through the shock wave is deflected sideways 
towards two polar and an equatorial accretion column, yielding a 
quadrupolar structure of the flow and the magnetic field, in con- 
trast to models with weaker fields, where the flow structures are 
dominated by higher multipolarity. 

Accretion proceeds along a few directions, viz. in these 
columns, at trans- and superalfvenic velocities, whereas the 
gas moves subalfvenically between these accretion flows (the 
blue regions in the right half of Fig.[T2li. In these regions, the 
Alfven points lie far out in the accretion/SASl region at more 
than 150 kilometres, whereas they are located roughly at the de- 
celeration radius in the downflows, i.e., at radii around 100 kilo- 
metres (see the borders between blue and red regions in the right 
halves of Fig.[T2li. While these regions exhibit some variations 
of their shape with time, they are present during the entire post- 
bounce time. It is also worth mentioning that is less than 
unity in most of the post-shock regions, meaning that perturba- 
tions of the magnetic field can indeed be described in terms of 
Alfven waves, whereas Ay varies between values greater and less 
than unity, i.e., the waves are either trapped by the flow or can 
escape. 
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Fig. 14. Comparison of the entropy per baryon of models 
sis -Bm (left half) and sl5-B12 (right) at t ^ 170 ms. 



Due to the stable quadrupolar structure of the flow, the dy- 
namics of an accretion column, as outlined in Sect. 14. 2.2] is im- 
portant for the evolution of this model. Magnetic field is accreted 
in the downflows, while it can rise against the accretion in the 
broad regions of slowly moving gas between. Similarly to model 
slS-Bll . 5, waves travel upwards along the polar axis from the 
PNS convection zone. The accretion flow is decelerated abruptly 
at a discontinuity at the lower end of the accretion column, and 
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Fig. 15. Snapshot of model sl5-B12 at t k 300 ms. The left and 
right halves shows the logarithm of the magnetic field strength 
and the entropy per baryon, respectively. 



an Alfven point is located there. At this point, we find a pro- 
nounced enhancement of the magnetic field. 

Fig. [13] displays a comparison of average profiles of the 
magnetic fields for models models sl5-B12, sl5-Bll, and 
slS-BlO, corrected for the scaling of the initial field strength. 
While the models with weaker initial fields (green and orange 
lines) show a distinct increase of the field strength at the outer 
boundary of the PNS convection zone, i.e., at r » 50 km, a simi- 
lar feature is present for the model with stronger field at a much 
larger radius of r ~ 100 km; in this case, the profile is described 
roughly by a uniform power law out to this point. The maximum 
field strength of model sl5-B12 is of the order of 10'^ G at the 
base of the convection zone. 

Between the accretion columns, strong magnetic flux sheets 
are located at the lateral and upper boundaries of the upwards 
flowing plumes. This is visible, e.g., in the large arcs of field 
connecting the accretion columns immediately behind the shock 
(see Fig.fTSI). Gas and magnetic field accreted in the accretion 
columns rises in these regions from the deceleration layer. On the 
other hand, rising magnetic flux can, upon reaching the bound- 
aries of the plumes, be entrained in the accretion flow and thus 
increase the accreted magnetic energy. 

In Fig. [141 we compare the entropy distribution of models 
sl5-Bl(9 {left half of the figure) and sl5-B12 {right half) at 
t * 170 ms. The stable downflows of the latter model transport 
matter of low entropy towards the PNS where it is gains energy 
by neutrino heating. In this model, heating leads to a higher en- 
tropy of the gas rising between the accretion columns than in 
model S15-B10. This is caused by a magnetic confinement of 
the bubble at small radii and thus a longer exposure of the gas to 
strong neutrino heating. 

The shock expands gradually until t » 340 ms when it 
reaches an outermost location of about r » 600 km. We show 
a snapshot of the model shortly before maximum shock expan- 
sion in Fig.[T5]and a space-time diagram of the model in Fig.[T6l 
At maximum expansion, the quadrupolar structure is replaced by 
more complex flows. The model exhibits a large-scale pattern of 
low-order multipoles. Then, expansion ceases and the shock re- 



treats again to r ~ 300 km, but it starts to expand again after 
t ^ 420 ms. At the end of the simulation, t x 500 ms, its max- 
imum radius is 600 km again. While large fluctuations of the 
shock radius can be found in other models as well, this model is 
characterised by a clearly stronger expansion. The shock shows 
appreciable deformations; however, we do not find a very large 
global asphericity. The axis ratio is around unity during the late 
phases of the simulation. This can be seen, e.g., in the space- 
time diagram in Fig. [161 the regions of essentially high values 
of the angular averaged entropy and magnetic field extend out 
to the maximum shock radius. In a model with a strongly de- 
formed shock, the maximum shock radius can be considerably 
larger than the radius out to which these variables assume high 
values. 

4.4. Comparison of global quantities of our models 

After the discussion of the magnetic-field amplification and of 
the structural differences of the simulated models in the previous 
sections, we conclude the presentation of our results by a com- 
parison of important global quantities characterising the overall 
evolution of our models. 

We evaluate the total mass and the kinetic and magnetic en- 
ergies contained in the gain region, i.e., in the region between 
the shock wave and the stable layer surrounding the PNS. In this 
region, the matter gains energy by reactions with neutrinos dif- 
fusing out of the PNS, and hot-bubble convection and the SASI 
operate. 

The mass in the gain region, mgain (shown for the non- 
magnetic model sl5-B0, one weak-field model (sl5-Bll), and 
the strong-field model sl5-B12 in the left panel of Fig.fTTI). in- 
creases during the first 50 milliseconds post-bounce and peaks at 
about 0.055Mq for all models. Afterwards, the mass decreases 
gradually, falling below O.O2M0 after f ~ 320 ms, in all but the 
strongest magnetised model. In the case of model s 1 5 -B 12 mgain 
instead levels off at a value of mgain ~ 0.04Mq. This different be- 
haviour reflects the different evolution of the shock wave. While 
the outermost shock position can reach a radius of ~ 500 km 
along the poles, the average shock radius oscillates around a 
value of Tsiiock S 200 km in the non-magnetised and the weak- 
field models. On the other hand, the shock expands to radii of 
300 - 500 km at all latitudes in the strong-field model. Hence, 
a higher mass can be affected by neutrino heating in the model 
with the strongest initial field investigated. Although our models 
do not allow us to draw any firm conclusions on the long-term 
evolution of the models, we may estimate tentatively from the 
evolution of mgain that model sl5-B12 is closest to a successful 
explosion. Further indication of a different behaviour of weak- 
field and strong-field models is provided by the total energy of 
the matter in the gain region, i.e., the sum of internal, kinetic, 
magnetic, and gravitational energies. We find that this total en- 
ergy of model sl5-B12 exceeds that of model sl5-BIS by about 
5 X 10^*' erg around f = 440 ms. 

The lateral kinetic energy in the gain region, e^j^.^^j^ (shown 
in the middle panel of Fig.[T71for the same models), indicates 
the strength of hot-bubble convection and of the SASI. In agree- 
ment with the evolution of mgain, the lateral kinetic energy of 
model S15-B12 exceeds that of the other models at late times 
(f > 300 ms) by a factor of roughly 2. Apart from this difference, 
we find a higher value of ef. . in model sl5-B 12 than in the 

^ Kin;gain 

Other models even at early times, f < 120 ms, when all models 
have a very similar mass in their gain regions. Thus, the lateral 
kinetic energy per unit mass, or the average lateral velocity, is 
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Fig. 16. Space-time diagram of model sl5-B12. The top and bottom panels show colour-coded the entropy per baryon and the 
logarithm of the magnetic field strength, respectively. Additionally, mass shells and density contours are plotted, using the same 
spacing as in Fig.|2] 



largest in model sl5-B12. This is a consequence of the peculiar 
quadrupolar flow pattern developing solely in this model. 

We compare the amount of field amplification in the gain re- 
gions of all magnetised models in the right panel of Fig.[T7] For 
all models, an initial rise of the total magnetic energy in the gain 
layer emag;gain is followed by a phase of a relatively constant total 
field energies. When taking into account the scaling of the ini- 
tial fields (see the dashed lines), we find a very similar amount 
of field amplification for all models. This result suggests that 
the amplification is dominated by kinematic processes and that 
the back-reaction of the field onto the flow is weak. We do find, 
however, that the average magnetic energy of model sl5-B12is 
systematically lower than expected from the scaling of the ini- 
tial fields (cf. Fig. [17] where the black dashed line is below the 



other Unes). This somewhat weaker amplification is consistent 
with the presence of large subalfvenic regions where the field 
is sufficiently strong to react back onto the flow. The Alfven- 
wave amplification enhancing the field in this model compared 
to models with weaker initial fields affects regions that are too 
small to shine up in this volume-integrated energy. 

Finally, we discuss the distribution of the magnetic energy 
in the different regions of the core as a function of initial field 
strength. We refer to Tab.|2] for a synopsis of our results. We 
compare the root-mean-square averages of the magnetic field 
strength in four domains: 

b\ii is the r.m.s. field in the PNS at densities above 10'^ gem"-'; 
bcm is the rm.s. field in the PNS convection zone; 
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Fig. 17. Left panel: Mass in the gain region of models sl5-B0 (orange line), sl5-Bll (green line), and sl5-B12 (black line) as a 
function of time. Middle panel: Lateral kinetic velocity in the gain region of models sl5-B6) (orange line), sl5-Bll (green line), 
and S15-B12 (black line) as a function of time. Right panel: Logarithm of the magnetic energy in the gain region of models sl5-B10 
(red line), sl5-Bll (green line), sl5-Bll . 5 (blue line), and sl5-B12 (black line) as a function of time. The solid lines show the 
total magnetic energy of the models and the dashed lines (for the models with Zjq > 10'° G) the same quantity scaled according to 
the initial strength, respectively. 



b^fo is the r.m.s. field in the hydrodynamicaUy stable region sur- 
rounding the PNS; 
^gain is the r.m.s. field in the gain (SASI) layer. 

Additionally, we compare the ratios of magnetic to internal or 
kinetic energy, yS'^'', in the same regions. 

For all domains, the mean fields increase monotonically with 
the initial field strength, b(). The maximum field is reached in the 
PNS. Since this region is hydrodynamicaUy stable during the 
simulated evolution, amplification is dominated by radial com- 
pression. From a weak initial field, typical pulsar fields of the 
order of 10'^ G can be achieved. With a more extreme ini- 
tial field of b() ~ 10'^ G in the pre-collapse stellar core, we get 
magnetar-like values of 10'^ G. Compared to the thermal energy, 
the magnetic field is negligible even for this case; it can, how- 
ever, dominate over the velocity field. Since there is no feedback 
of the field, the final fields in the PNS core scale very well with 
bo. 

In the PNS convection zone, the average field is somewhat 
weaker than in the core of the PNS. As discussed above, the 
convective velocities can become subalfvenic for the strongest 
fields, corresponding to a magnetic energy exceeding the kinetic 
energy, i.e., > I; this is the case for model sl5-B12. In such 
a situation, the efficiency of the field amplification is limited. 
Consequently, the scaling of the r.m.s. field with bo breaks down 
for ^0^3x10" G. The field found for model sl5-B12, /jcnv ~ 
4.7 X 10^"^ G, should be close to the maximum field that can be 
achieved by the PNS convection excited in this core. 

Similar results can be found in the hydrodynamicaUy stable 
layer between the PNS convection zone and the gain region. The 
r.m.s. magnetic field strength is weaker by a factor of ~ 3. Again, 
we find a ratio of magnetic to kinetic energy in excess of unity 
for the strongest initial field considered and a deviation of the 
final magnetic field strength from a simple scaling with the initial 
field strength. 

The r.m.s. field strength in the gain layer is between 3.6 x 
10" G for model sl5-BlQ and 1.2 x 10^^ G for model 
S15-B12. In all models, is below unity, though. As discussed 
in Sect. l4.2.2l and Sect. 14. 31 trans- and subalfvenic regions are 
present and have an important influence on the dynamics and 



the field amplification. Because the magnetic field strength de- 
creases more slowly with radius than the internal energy, /3' can 
be as large as one per cent here. 

Finally, we note that for all models by far most (of the or- 
der of two thirds) of the magnetic energy is contained in the 
PNS convection region. The magnetic energies in the stable layer 
are of the same order of magnitude as those in the gain region, 
though the exact value of bg^in is variable and depends strongly 
on the size of the gain region. 



5. Summary and conclusions 

We have studied (some of) the processes leading to the amplifi- 
cation of the magnetic field in a non-rotating stellar core during 
the collapse and post-bounce accretion phases of a supernova. In 
non-rotating stars, a variety of amplification mechanisms con- 
sidered in previous works, e.g., winding of the field by differen- 
tial rotation or the MRl, are not viable. Instead, convection and 
the SASI may constitute (small-scale) dynamos. Additionally, 
Alfven waves travelling against the gas flow may be ampli- 
fied once they reach an Alfven point where the (co-moving) 
Alfven velocity equals the gas velocity, a condition that could 
be fulfilled in the accretion flow onto the PNS. If efficient am- 
plification occurs, the field may be able to affect the dynamics 
of the core, e.g., by altering the geometry of SASI flows or by 
energy dissipation through Alfven waves in the upper layers of 
the hot-bubble region. 

Reducing the complexity of the problem, simplified mod- 
els, e.g., one-dimensional simulations with assumptions for 
the excitation, pro pagation and dissipation of Alfven waves 
(Suzuki et al."2008'), toy models of Alfven waves in decelerat- 
ing flows (Guilet et al. 2010), and 2D and 3D MHD simulations 
without neutrino transport (lEndeve et al.l 12010 ) have demon- 
strated that these effects could, in principle, be relevant. On the 
other hand, to study their evolution under less idealised condi- 
tions, in particular their interplay with a highly dynamical back- 
ground, more calculations of self-consistent models are required, 
i.e. multi-dimensional MHD including a treatment of neutrino 
transfer through the stellar core. 
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To this end, we have performed axisymmetric simulations 
of the collapse and the post-bounce evolution of the core of a 
star of 15 solar masses possessing a purely poloidal initial field. 
Using a new code for neutrino-magnetohydrodynamics, we have 
solved the MHD equations coupled to the system of two-moment 
equations for the neutrino transport; the closure for the moment 
equations was provided by an analytic variable Eddington fac- 
tor. We included descriptions for the most important reactions 
between electron neutrinos and antineutrinos and the stellar mat- 
ter, viz. nucleonic and nuclear emission and absorption and scat- 
tering off nucleons and nuclei. Compared to cuiTent state-of-the- 
art techniques such as Boltzmann codes, our approach ignores 
muon and tau neutrinos and is less accurate in treating neutrino- 
matter interactions, but needs much less computational resources 
and possesses full two-dimensionality including velocity effects. 
We are able to reproduce the basic features found in detailed sim- 
ulations, e.g., the stagnation of the prompt shock wave, PNS and 
hot-bubble convection, and the S ASI activity. Hence, our models 
allow for a fairly reliable assessment of the main MHD effects in 
self-consistent SN core models; we discuss the remaining major 
limitations below. 

The principal results of our simulations and the main conclu- 
sions to be drawn from our results are: 

1 . The magnetic field is amplified kinematically by the turbu- 
lent flows developing due to convection and the SASI. The 
amplification factor does not depend on the initial field un- 
less it reaches, starting from a high initial value, equipartition 
with the kinetic energy, the maximum possible energy attain- 
able. The maximum field we observe is a few times 10''' G at 
the base of the PNS convection layer. In the PNS convection 
zone, the radial profile of the field strength can be approxi- 
mated by a power law, \b\ oc r^^. 

2. When falling through the non-spherical supernova shock, the 
magnetic field is bent by lateral flows, creating a component 
parallel to the shock. These perturbations are advected to- 
wards the PNS convection zone. The field accumulates there, 
leading to a layer of strong magnetic fields. Since the pertur- 
bations are associated with a strong lateral component of the 
field, this layer is dominated by the 0-component of the field. 
Hence, Alfven waves propagate at constant radius rather than 
upwards. Therefore, we do not find an Alfven surface in 
the accretion flow, although there are sub- as well as super- 
alfvenic regions in this layer. This limi ts the efficiency of the 
amplification mechanism proposed bv lGuilet" 

etaDdmS- 

3. Conditions for the latter effect are more favourable if accre- 
tion occurs through a rather stationary column. In our ax- 
isymmetric models, this is the case mostly along the polar 
axis where the geometry enforces radial fields and flows. An 
Alfven surface forms in the radial field of the accretion col- 
umn at a radius depending on the initial field strength. For 
fields of b() > 10" G initially, we observe that perturbations 
created in the PNS convection zone propagate along the field 
upwards into the accretion column. Interacting with the ac- 
cretion flow, they lead to a steepening of the profile of the 
accretion velocity near the lower end of the accretion colum 
to a discontinuity where a strong enhancement of the mag- 
netic field takes place. 

4. The strongest initial field we have investigated, bo - lO'^ G, 
is able to shape the post-shock flow. It favours the formation 
of a flow dominated by low-order multipole modes, viz. a 
quadrupolar pattern of accretion columns at the poles and 
near the equator. In the stable polar columns we observe, as 
in the case of lower fields (bo = 10" G), an interaction of 



perturbations travelling upwards from the PNS convection 
zone with the accretion flow. This gives rise to the forma- 
tion of a layer of very strong field. Also in the convective 
and SASI layers behind the shock, the field strength grows 
with the magnitude of the initial field. In this model with its 
particular flow topology, we find the most pronounced shock 
expansion of all models. 

In summary, our results suggest that magnetic field amplifi- 
cation to interesting strengths can efficiently take place during 
the stellar core collapse even in the absence of rotation. In ad- 
dition to an enhancement due to compression by the radial col- 
lapse, we find that non-radial fluid flows associated with convec- 
tion and SASI activity, and the interaction of Alfven waves in the 
accretion flows can amplify the initial iron-core fields. Starting 
with 10^-10'° G as pr edicted by present stellar evolution mod- 
els dHeger et aLlliOOSh . field strengths of typical pulsars (10 - 
10'^ G) can be reached. Magnetar fields of lO'^'-lO'^G result 
when the progenitor core possesses a pre-coUapse field between 
a few 10" G and 10'^ G. Only in the latter case, the fields around 
the nascent neutron star obtain dynamical importance and might 
have an influence on the supernova explosion mechanism. 

Though offering some insight into the magnetic-field evo- 
lution in non-rotating magnetised cores, our study has several 
important limitations: 

1 . We have used a new treatment of neutrino transport and sim- 
plified neutrino-matter interactions constraining ourselves to 
electron neutrinos and antineutrinos. While the accuracy of 
the transport compared to other schemes has to be assessed in 
a separate study, and our set of neutrino reactions is certainly 
insufficient for high-precision supernova physics, we do not 
deem this a problem for the presented investigation because 
we are able to capture the most important dynamical effects 
in a supernova core. 

2. To save computational costs, we have restricted ourselves to 
axisymmetric simulations. In the light of the anti-dynamo 
theorems, this is a severe limitation leading to a wrong, 
possibly too low, level of amplification of the field in the 
turbulent regions. Furthermore, our models do not allow 
for the development of shear layers associated with non- 
axisymmetric spiral modes of the SASI, which may also be 
a site of effici ent amplification of the magnetic field (see also 
lEndeveet"ani2() 10). 

3. Moreover, axisymmetry restricts the dynamics of the accre- 
tion flows, favouring the development of very stable accre- 
tion columns along the poles. As we have seen, field am- 
plification shows very distinct features in and below these 
columns. We presume that the dynamics in three dimensions 
is dominated by what we have seen in off-axis accretion 
flows, i.e., less stable Alfven surfaces and less efficient am- 
plification of Alfven waves. Very strong initial fields of the 
order of bo = 10'^ G lead to high field strengths that are able 
to dominate the post-shock accretion flow. This may estab- 
lish coherent, stable accretion columns. For such fields, our 
axisymmetric results may hence be a decent approximation. 

4. Turbulent field amplification may be very sensitive to dis- 
sipation coefficients, physical and numerical. Our models, 
based on ideal MHD, neglect dissipation by physical vis- 
cosity and resistivity, but are computed on relatively coarse 
numerical grids corresponding to excessive numerical dis- 
sipation. Therefore, we are not able to follow the turbu- 
lent (inverse) cascades of magnetic and kinetic energy and 
helicity covering many orders of magnitude in wave num- 
ber in a supernova core. The effect of insufficient reso- 
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lution on Alfven waves is probably less serious though 
their wave number should increase as they approach the an 
Alfven point, requiring enhanced resolution. Simulations of 
cores at a resolution corresponding to numerical viscosity 
and resistivity below the physical ones are by far too expen- 
sive today, and will remain so for a long time. To tackle this 
difficulty, a combination of different approaches would be 
desirable, viz. global direct numerical simulations of the core 
with drastically enhanced physical transport coefficients, and 
simple sub-grid models for MHD turbulence based on ide- 
alised local simulations neglecting most aspects of, e.g., neu- 
trino physics. We are, however, aware of the lack of reli- 
able sub-grid models for MHD at present, obstructing further 
progress in this direction. 

Apart from these methodological shortcomings, open phys- 
ical questions are, e.g., the influence of the progenitor on the 
establishment of certain patterns in the accretion flow, effects of 
very strong initial fields, and the influence of slow rotation of the 
core on our findings. In particular the last issue may prove inter- 
esting as it would enable a large-scale dynamo. We defer these 
questions as well as the more technical problems listed above to 
further investigations. 
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Appendix A: Neutrino transport 

Most results on magnetic effects in supernovae draw from more 
or less idealised investigations, e.g., analytic calculations or sim- 
plified models focusing on one particular aspect of the problem 
while neglecting some of the complex physics of the explosion 
and, thus, limiting their range of application. Due to the enor- 
mous computational requirements of an accurate treatment of the 
transport of neutrinos and their interaction with matter, simplifi- 
cations are most commonly made in this sector Not discussing 
local simulations focusing on a high-resolution modelling of a 
small part of the flows, we can distinguish between different ap- 
proximations: 



1. Some studies (e.g.. lObergauUngeret al.ll2006HMikami et al.l 
l2008h neglect the effect of neutrinos completely and use only 
a simplified, e.g., polytropic, equation of state. Typically, 
these models exhibit prompt explosions, limiting their valid- 
ity to the immediate (few milliseconds) post-bounce phase. 
Even without any treatment of neutrino s, it is possible to sim - 
ulate the development of the SASI (e.g.. lBlondin et al.'2003). 

2. A way to model the dynamics qualitatively correctly while 
avoiding the complications of a full neutrino-transport 
scheme is the use of simple heating/cooling or deleptoni- 
sation source terms. Such schemes prescribe the loss and 
gain of energy or electron number by cooling and heat- 
ing terms that depend on the local thermodynamic state 
of the matter (e.g., density, temperature, electron fraction). 
Depending on the details of the model, the can reproduce, 
e.g., the correct deleptonisation behaviour during collapse 
(e.g., Liebendorfer 2005) or t he development o f the SASI 
(iFernandez & Thompson>2009HFernandezl20I0l) . 



A leakage scheme accounts for the optical depth of the re- 
gion emitting neutrinos is taken into account in the cooling 
terms. The results of such simulations can be s imilar to that 
of models with purely local coohng terms (e.g.. lKotake et al.l 
i2004 . 

Despite the large computational effort, various methods for 
energy-dependent neutrino hydrodynamics have been used; 
state-of-the-art codes include, e.g., the isotropic diffusion 
source approximation, a scheme decomposing the neutrinos 
into a trapped, diffusing, and a freely streaming component 
dLiebendMer et al.ll2009l) . flux-Hmited diffusion (FLD), or 
solve the Boltzmann equation for the neutrino phase-space 
distribution. Among the codes based on the FLD method, 
there is a wide variety of approaches to the treatment of 
mu lti-dimensio n al effe c ts and the coup li ng of energy bins: 

- iLivne et all (|2004 : iBurrows etall (l2()07h use a multi- 
group FLD code without energy-bin coupling and no ve- 
locity terms i n the m o ments equa ti ons, 

- IWalder et all (l2005h : lOtt et al.l (l2008h use a two- 
dimensional multi-angle, multi-group FLD code without 
energy-bin coupling and no velocity terms in the moment 
equations, 

- the simulat i ons of iMezzacappa et al.l (l2007h : 
lYakunin et al.l (1201 Oh employ a ray-by-ray-plus multi- 
group FLD method that includes the energy-bin coupling 
terms and the velocity terms up to order (9(v/c). 

Similar distinctions can be made for the codes based on 
Boltzmann solvers: 

- Mezzacappa & Bru eniil (Il993h; IMezzacappa et al.l 
(120011) : iLiebendorf er et all (l2004 l) solve the one- 
dimensional Bol t zmann equation by an 5 „ technique, 

- iRampp & Jankal (l2002h employ a one-dimensional two- 
moment-closure scheme with a closure relation obtained 
from the solution of a model Boltzmann equation, 

- this scheme wa s generalise d to a ray-by-ray-plus ver- 
sion by Buras et al.l (l2006h . accounting for all terms 
coupling energy bi ns and all velocity-dependent terms. 
iMiiller et al.l (l20I0h developed a 2D general relativistic 
hydro sol ver with a r ay-by-ray-plus neutrino transport 
similar to lBuras et al .' (2006). 

These schemes provide the best models of the complex flow 
patterns of a supernova and allow for a reliable distinc- 
tion between successful and failed supernovae and an ac- 
curate prediction of the explosion energy. Testing the lim- 
its of currently available supercomputers, two-dimensional 
axisymmetric simulations are possible, but the step to three- 
dimensional models may require a profound change of al- 
gorithms and programming models to take full advantage of 
the next generation of computers. 



Neutrino transfer describes the neutrinos in terms of their 
phase-space distribution function, f(t, x, e, n), a function of time, 
f, position, X, neutrino energy, e, and propagation direction, n 
{\n\ = 1). The dimensionality of this function is the main rea- 
son for the high computational requirements of RT. It can be re- 
duced by evolving the set of neutrino moments. My - J dO n*/, 
k = 0, oo, rather than the distribution function; the first mo- 
ments have an immediate physical meaning: k = 0, 1,2 corre- 
spond to the neutrino energy density, momentum density, and 
pressure tensor, respectively. 

In principle, the full infinite series of moments is required 
for an accurate representation of /. In practise, the series is trun- 
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cated at some level, ^max, retaining the conservation laws for the 
first few moments. Since the equation for the A:* moment. 



(A.l) 



involves at least the next moment (order A: + 1), a closure is re- 
quired to relate the first moment beyond truncation, mJ°'"^', to 
the moments of order 0, ...,fcmax (the source term, Q\„ accounts 
for all interaction terms and for velocity-dependent terms such 
as advection, compression, or Doppler shifts). A common choice 
is the 0*-moment system, the equation for the neutrino energy, 
with the first moment given by the (flux-limited) diffusion flux. 
We go one step beyond this and evolve the P'-moment system. 



d,Ey + VFy 
d,Fy + VPy 



Ql 



(A.l) 
(A3) 



closing the system by assuming a simple analytic form of the 
2"'' moment, Py = Py(Ey, Fy). In a one-dimensional system, one 
scalar function is su fficient to define th e closure, the Eddington 
factor Py. Following lAudit et al.l (l2002h . we construct a tensorial 
closure for multi-dimensional systems from the one-dimensional 
Eddington factor: 



p'J = 



5" + 



3py - 1 FlF; 
2 F,2 



(A.4) 



This expression combines a diagonal part, accounting for 
isotropic diffusion, with a contribution accounting for free 
streaming in the direction of the neutrino flux. In the diffusion 
limit, the Eddington factor is py - I /3, and the pressure tensor is 
given by the isotropic part alone, P'J = Ey/36'^, and in the free- 
streaming limit, Py = I, only the second part contributes, the 
total neutrino pressure is Ey and P'J is peake d in the propagation 
direction of the neutrinos. Pons et al. '(2000) have demonstrated 
the applicability of standard HRSC methods to this (hyperbolic) 
system and compared different suggestions fo r the closure. In 
our simulations, we use the Eddington factor bv lMinerbol([T978h 
To avoid the evaluation of the Langevin function in the expres- 
sion for the closure, we apply a simple polynomial approxima- 
tion given by Cernohorsky & Bludman (.1994.) : 



6/;), 



(A.5) 



where fy - Fy/cEy is the flux factor of the neutrinos, i.e., 
the neutrino flux normalised to the maximum flux allowed by 
causality. 

Exploiting the hyperbolic character of the RT system, we 
use HRSC methods for the moments equations, viz. the same 
high-order reconstruction as in the MHD subsystem, either the 
Lax-Friedrichs or the HLL approximate Riemann solver, and an 
explicit Runge-Kutta time integrator. This approach meets its 
limitations when the Peclet number, Pe - x^^^^ i-^., the opti- 
cal depth of one grid cell of width 6x (c is the speed of light), 
exceeds one. Then, the interaction source term of the neutrino 
fluxes, —cx^Fy, becomes stiff" (the characteristic time scale being 
~ 1/(0^)); enforcing a time step much below that of the hyper- 
bolic MHD or transport parts, (5f)hyp ^ ^x/cm ax^ for ^ given max- 
imu m characteristic sp eed c^ax. Similarly to iPons et al.l (l2000l) 
and lAudit et al.l (|2002|) . we modify the Riemann fluxes of the O* 
moments to reproduce the correct parabolic difiiision limit and 
use an implicit time integrator for the interaction source term 
of the first moment, thereby avoiding a time-step restriction be- 
low (5f)hyp. A stiff part of the system on its own if cx^^6t > 1, the 



neutrino-matter interaction source terms in the O^-order moment 
equation are treated implicitly as well. 

Including contributions due to the fiuid velocity of or- 
der 0(v/c), we solve the equations for the neutrino moments 
in the co-moving frame of reference, summarised, e.g., by 

Im 

unier & Weaven is the neutrino energy): 



Jk 



d,Ey + V ■ (EyV + Fy) - LoVjVkd^P^ 

d,F'y + Vj(Fy +P'J) + F'yVjvj = Si 



(A.6) 
(A.7) 



In these equations, we have separated the total source terms, 
Qy'^ (Eq. ( IA.2I 1 and Eq. ( IA.3l l), into source terms due to neutrino- 
matter reactions, 5v'', and velocity-dependent terms. We ensure 
the simultaneous conservation of neutrino energy and number 
density by evolving a duplicate set of neutrino moments for both 
conserved variables. There are alternatives to this approach, e.g., 
using a consistent d iscretisation of the neutrinos in energy space 
(lMulleretal.1120101) . 

Striving for simplification of our modelling, we include a 
reduced set of neutrino-matter interactions only. It consists of 

1 . emission and absorption of electron neutrinos by neutrons, 

2. emission and absorption of electron anti-neutrinos by pro- 
tons, 

3. elastic scattering of all neutrino fiavours off nucleons, 

4. emission and absorption of electron neutrinos by heavy nu- 
clei, 

5. coherent elastic scattering of all neutrino fiavours off heavy 
nuclei. 

The reaction ra t es and opacities are implemented following 
Rampp & Jankal (|2002). Neglecting pair processes, we do not 
create jj and r neutrinos. Furthermore, we neglect inelastic scat- 
tering of neutrinos off leptons. From these reactions, we can con- 
struct the source terms of the 0* and P' moments. 



S^y = cx"(Ey'-Ey) 
Si = -cx'Fy, 



(A.8) 
(A.9) 



where El'^ is the equilibrium (Fermi-Dirac) value of the 0* mo- 
ment. The opacities for the 0* and moments include absorp- 
tion and emission processes only and absorption, emission, and 
scattering processes, respectively. 

The number of neutrinos in a certain energy bin is lim- 
ited by the Pauli exclusion principle. If a bin is populated by 
a number of neutrinos close to this limit, interactions, in par- 
ticular inelastic scattering, redistributing the particles in energy 
space, are getting phase-space blocked such that a violation of 
the Pauli principle is avoided. Numerical schemes solving the in- 
teraction terms and the spatial derivatives of the transport equa- 
tions time-implicitly, obey this constraint naturally. Since we are 
using an explicit time integrator for the spatial derivatives in- 
stead, we have to enfo rce this constraint manually similarly to 
iSwestv& Mvral (120051) . 
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